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1 • Justification and Objectives 

Among the tools used in passive remote sensing of Earth resources in the visi- 
ble and near-infrared spectral regions are measurements of spectral signature and 
k^i tect ional reflectance functions (BDRFs). Determination of surface properties 
using these observables is complicated by a number of factors, including (a) mixing 
of surface components, such as soil and vegetation, (b) multiple reflections of 
radiation due to complex geometry, such as in crop canopies, and (c) atmospheric 
effects. Differences in surface conditions, spectral coverage, viewing and illumi- 
nation geometry, and atmospheric conditions make comparison of observations obtained 
in the laboratory, in the field, and from aircraft and spacecraft extremely diffi- 
cult. In order to bridge the diversity in these different approaches, there is a 
need for a fundamental physical understanding of the influence of the various effects 
and a quantitative measure of their relative importance. In particular, we consider 
scene complexity effects using the example of reflection by vegetative surfaces. 
The Interaction of sunlight with a crop canopy and interpretation of the spectral 
and angular dependence of the emergent radiation is basically a multidimensional 
radiative transfer problem. In the first st?.p of the radiation historv, sunlight 
is directly and diffusely transmitted downward through the Earth’s atmosphere. The 
downward radiation field then interacts with the plant canopy. Each element of the 
canopy (leaves, stems) reflects (and transmits) light according tc Its Intrinsic 
bidirectional reflectance (and transmittance) properties. The complex canopy geo- 
metry, underlying soil cover, and presence of diffuse as well as collimated illumi- 
nation will modify the reflectance characteristics of the canopy relative to those 
of the individual elements. Finally, upwelling radiation emergent from the canopy 
will t>c directly and diffusely transmitted \>y the atmosphere to space, .and .some 
contamination of the* spectral signature with radiation reflected from neighboring 
surface regions will occur due to scattering associated with the adjacency effect. 
The intent of this research program is to develop some of the tools needed In order 
to achieve a quantitative understanding of such phenomena, using a combination of 
experimental and theoretical methods. 

2* Technical Approach 

The tools which are needed in order to study this radiative transfer problem 
include (a) knowledge of the bidirectional scattering properties of Individual 
canopy elements, primarily leaves, as well as of the canopy as a whole, (b) a 
model describing the interaction of radiation with the canopy, and (c) a radiative 
transfer algorithm for describing the interaction of a nonuniform radiation field 
with the atmosphere. In order to obtain realistic leaf bidirectional reflectance 
functions for incorporation into canopy and scene models, we are using laboratorv 
spect rogoniometry to atudy the intrinsic scattering properties of Individual leaves. 
These measurements provide a data set which obviates the need for nonphysical model 
assumptions, such as Lambertian leaf BDRFs. We have been using a goniometer located 
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at the University of California, Los Angeles to obtain spectra of Individual leaves 
for a variety of Illumination and viewing geometries. The ln8t ^°\ 1 ? c e e r d a “ t 3 ^ 
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f 3 r f fer !! d 5° 3 f Model K Figure 2 shows Che results of .alculations for 
Model 1 using the Fourier transform technique, along with selected points from the 
Monte Carlo results. The total upwelling Intensities are shown as the sum of the 
diffuse (multiply scattered) radiation field and the direct field, which consists 
of photons transmitted to space without being scattered. The agreement between 
the two methods appears to be excellent. We also show at large distances from the 
boundary the Intensities calculated using the one-dimensional matrix operator 
method for uniform albedos of 0.0 and 0.6, demonstrating that the Fourier transform 
results approach the expected values away from the boundary. In the vicinity of 
the boundary, the adjacency effect Is apparent. Due to scattering within the atmo- 
sphere, the upwelling Intensities on either side of boundary are affected by photons 
reflected from the surface on the opposite side of the boundary and then scattered 
nto the line of 3ight. In effect, the atmosphere "blurs" the boundary with a point 
spread function whose width is a few kilometers. The Fourier transform technique 
s useful in that its underlying mathematical formalism demonstrates how the atmo- 
sphere £ l}ters high spatial frequency Information from the upwelling radiation 
Held. This filtering can be described by the computation of atmospheric modulation 
transfer functions (MTFs). Figure 3 Is a comparison of MTFs calculated using our 
technique for a Rayleigh atmosphere and an aerosol laden (hazy) atmosphere, each 
of total optical depth 0.3 and 5 km scale height. These MTFs indicate that the 
more highly forward scattering haze phase function results in a greater observabll- 
ity of the adjacency effect. 


Another area in which the 3-D radiative transfer algorithm is proving useful 
Is In the development of a technique for determining atmospheric transmittance 
from surface Images acquired at off-nadlr view angles. The basis of the method Is 
as follows; A3 shown In Fig. 2, the total upwelling radiance is the sum of a direct 
and diffuse component. Radiation directly transmitted to space Is proportional to 
the surface albedo distribution, and a sharp discontinuity in the surface albedo 
gives rise to a similar discontinuity in the direct field. The diffuse field, on 
the other hand, while exhibiting a qualitative similarity to large scale variations 
in surrace albedo, exhibits a much smoother character, owing to the filtering effects 
of the atmosphere, as discussed above. The direct field, which is also proportional 
to atmospheric transmittance exp(~r/p), where r is the optical depth and p is the 
cosine of the view angle, thus has a very different spatial character from the 
diffuse field. If we compute the spatial Fourier transform of upwelling radiances 
within an image, the spatial "smoothness" of the diffuse field means that it will 
.make-a .negligible -contributidn-at high spatial frequencies. Thus, if we eliminate 
xow frequency Fourier components we are left with a measurement that is directly 
proportional to the atmospheric transmission factor exp(- T /p). By obtaining images 
of the same surface region at several view angles, it should in principle be possible 
to solve for x, provided non-Lambert lan surface effects can be minimized. Our cal- 
culations suggest that acquisition of images in the vicinity of 60° off-nadir should 

be sufficient to accomplish this and recover atmospheric transmittance with an accu- 
racy of 5-10%. 

Significance of Results and Future Directions 

Many investigators have studied the spectral properties of a wide variety of 
crop samples. The spectral reflectance of the foliage is modified by changes in 
the internal structure of the leaves associated with various forms of stress, such 
as water or mineral stress. The specularity of the reflected radiation increases 
at wavelengths where increased absorption occurs, presumably due to reduced depth 
o. penetration of the incident radiation into the leaf resulting in fewer multiple 
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reflections between layers. Therefore, we expect changes in leaf structure result- 
ing from stress to affect the angular as well as spectral reflectance properties. 
Thus, it is important to incorporate realistic (non-Laabertian) leaf bidirectional 
reflectance functions into canopy and scene models. The spectrogonionsetric leaf 
measurements will provide such data* 

The Fourier transform radiative transfer algorithm which we have developed is 
a useful tool for studying the effect of the atmosphere on high spatial resolution 
imaging of the Earth's surface. As discussed above, the sharpness of the adjacency 
effect in the vicinity of albedo boundaries is dependent upon the vertical distribu- 
tion of the atmospheric scatterers. We are using our radiative transfer method to 
perform a systematic study of the variation of atmospheric point spread functions 
with aerosol vertical distribution, phase function, sin v^e scattering albedo, and 
optical thickness. The results of such calculations will permit evaluation of the 
possibility of observing atmospheric blur directly in high resolution images of 
the Earth's surface, with the goal of determining aerosol properties. Conversely, 
these calculations will lead to techniques for removing spectral contamination 
associated with the adjacency effect. 

Further work is necessary in order to demonstrate the feasibility of implement- 
ing the multiple view angle optical depth retrieval technique on a space platform. 
The major uncertainty at present is the manner in which the unknown angular reflec- 
tance properties of the surface should be taken into account. For the model BDRF we 
used in our simulations, off-nadir imaging near 60° appears appropriate to compensate 
for lack of a priori knowledge of the shape of the surface BDRF. Whether this range 
Is appropriate for natural surfaces can only be determined by experimentation in the 
field. A number of portable field spectrometers are available for obtaining BDRF 
measurements. The success of this technique would be significant because it provides 
a means of monitoring tropospheric opacity, which is important not only for atmo- 
spheric studies but also because it is one of the most important parameters govern- 
ing the effect of the atmosphere on surface images. An extension of the optical 
depth retrieval technique is to near infrared spectral regions where water vapor 
absorption occurs. In Its current form, the technique relies on the assumption that 
atmospheric transmittance variesias exp(-T/y). In the water vapor bands, any finite 
bandwidth spectral channel will incorporate many individual lines and this law may 
not be followed. In order .to study the view angle dependence of transmittance in 
* ’ such* regions, *we‘ are ‘developing aradiative transfer algorithm which- permits * calcu- 
lation of broadband radiances for channels containing many absorption lines in tne 
presence of aerosol scattering. 

Finally, we plan to Investigate the applicability of our 3-D radiative trans- 
fer method to scattering in plant canopies. In principle, our algorithm is appli- 
cable to heterogeneous media and has the advantage of providing a true solution 
to the multidimensional radiative transfer equation. In addition, solutions are 
obtained at many view and illumination angles and azimuths, thus providing the 
full canopy BDRF. Variation of the BDRF of the individual leaves will permit 
study of the effect of changing reflectance properties of the Individual canopy 
elements on the canopy angular and spectral reflectance is a whole. Such calcula- 
tions will mesh nicely with the goal3 of the experimental portion of this investi 
gation. 
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A 1-D MATRIX OPERATOR 
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